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We calculate the ground-state energy and other physical properties of the hydrogen atom inside a 
spherical box with an impenetrable wall. We apply the variational method and perturbation theory 
f and compare both approximate results. We show that the total, kinetic and potential energies for 

the moving-nucleus model are greater than those for the case in which the nucleus is clamped at 
' the box center. 

^ ■ I. INTRODUCTION 

in 



Atoms and molecules confined into boxes of different shapes and permeability proved to be suitable simple models 
for the study of the effect of the environment or of high pressure. In the later case the pressure is given by p = —dE/dQ 



where E is the energy and ft the volume of the confining box. Those models enable us to simulate the effect of pressure 
on several properties of atoms and molecules, such as, for example^ their polarizability[l|. There is a vast literature 



on the subject and the reader may find suitable reviews elsewhere[ 

In all those cases the nucleus was clamped at some chosen point within the box. For example, in the case of an 
impenetrable spherical box we know that the energy is lowest when the nucleus is clamped at origin and increases 
as the nucleus approaches the surface 7}. This behavior comes entirely from the interaction between the electron and 
the wall and does not tell us anything about the effect of the hard surface on the nucleus. 

If we have an atom in a real environment we should assume that the nucleus also interacts with that environment. 
Therefore, it seems reasonable to add that interaction explicitly into the model. The simplest one is that in which 
both the nucleus and electron are affected exactly in the same way by the impenetrable spherical surface. In such 
a case we know that the nucleus cannot have zero kinetic energy because of its interaction with the surface and, 
therefore, the clamped-nucleus approach may not be the most realistic one 80. 
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It seems more reasonable to assume that a repulsive force on the electron should be an attractive one on the nucleus 
and conclude that it is more realistic to choose different boundary conditions for each particle at the box surface. 
However, as a first approximation to more elaborate models it seems sensible to start with an impenetrable wall for 
both particles and compare the results of this model with the clamped-nucleus approach. 



II. THE MODEL 



The Hamiltonian operator for a nonrclativistic hydrogen atom is 

H = f + V 

f = — — v 2 --^— V 2 

2m e e 2m„ 11 

V(r) = (1) 

where m e and m n are the masses of the electron and nucleus located at r e and r„ with charges — e and e, respectively, 
r = |r|, r = r e — r n , eo is the vacuum permittivity and V 2 denotes the Laplacian in the coordinates indicated by the 
subscript. 

If the atom is confined to a spherical box of radius R with an impenetrable wall we can formally write the Hamil- 
tonian operator in atomic units as 

H = f + V(r) + U(r e ,r n ) 

mh"'""' 8 P) 

oo if r e , r n > R 

where m n = 1836.15267261 and the eigenfunctions should vanish when either r e = R or r n = R. 



It is well-known that the virial theorem for such a system is given by 



dE „ dE I ~ \ 

where £1 is the volume of the spherical box. This expression gives us the pressure on the hydrogen atom p = —dE/dfl 
in terms of the kinetic and potential energies. When R — >• oo we recover the virial theorem for the free hydrogen atom 

2 (f r ) + (V) =0,f r = -tt^V 2 (4) 

where T r is the kinetic energy for the relative motion and m = m n /(m n + 1) is the reduced mass in atomic units. 



III. METHODS OF CALCULATION 
A. The variational method 



In order to calculate accurate eigenvalues and eigenfunctions of the hy drog en in a box we resort to the generalized 
Hylleraas basis set that proved useful for three-body Coulomb systems 12j4l5j • For example, for the S states we 
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choose a trial function of the form 

N 

V(r e ,r n ,r) = (l - ^) (l - ^) £ ^^^-"^.-^-7^ (5) 

fc=i 

with linear Cfc and nonlinear ak,fik,lk variational parameters. The explicitly correlated character of these functions 
ensures accurate energies for the ground and low excited states of free three-body atomic and molecular species with 
relatively few terms in the expansion. The effect of the Hamiltonian operator (|2|) on this variational function is simply 
given by 

H - 1 ( ^ I 2 9 1 r e ~ r n + 7-2 9 2 \ 

2 \dr 2 r e dr e r e r dr e dr J 

1 / d 2 2__8_ rl-r 2 e +r 2 d 2 \ 
2m„ ydrl r n dr n r n r dr n dr ) 

m n + 1 / d 2 | 2 d \ 1 

2m n \ dr 2 r dr J r 

A recent pedagogical approach to the confined hydrogen atom with a moving nucleus proposed the simple ansatz 

V [r e , r„) = N( 7 , R) (l - |) (l - ^) e"^ (7) 

where 7 is a variational parameter and iV(7, R) the appropriate normalization factorQ, Since this function yields 
the correct result when R — > 00 then it provides accurate ground-state energies for large and moderate values of R. 
In this paper we propose the straightforward generalization 

N 

Are, r„) = (l - T -£) (l - e -*.-^--rr ^ Cfc r^<V* (8) 

fe=i 

that may be suitable for the ground and excited 5* states. 

Notice that the cut-off function (1 — r e /R)(l — r n /R) tends to unity as R — !> 00 so that the resulting asatz © will 
be suitable for the free system. We also expect that under such condition a and j3 vanish as well as the coefficients 
Cfc of those terms with nk ^ and mk 0. 

For comparison we will also consider the hydrogen atom with the nucleus clamped at origin. A suitable trial function 

is 

JV 

\ P r Sr 

Rj 



<P CN (r) = (l-^)e- Sr Y,djr j (9) 
3=0 

where S and dj are nonlinear and linear variational parameters, respectively. This trial function yields the exact 
S states of the free hydrogen atom when R — > 00. For brevity we refer to the moving-nucleus case (MNC) and 
clamped-nucleus case (CNC) from now on. 



B. Perturbation theory 



In the strong-coupling regime R — > the Coulomb interaction is entirely dominated by the kinetic energy and the 



problem is almost separable j8-ll0l|. In the first approximation we have two particles moving independently within 
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the box and we can apply Rayleigh-Schrddinger perturbation theory by splitting the Hamiltonian operator into the 
zero-order or reference 

Hv = -\vl--^-Vl + U{r e ,r n ) (10) 

and perturbation H' = —l/r parts. 

The unperturbed eigenfunctions are given by 

€L'lw( r ei r J = Xnlm(r e )Xn'l'm'(r n ) (11) 

were 

Xnlm(r) = Rni(r)Y lm (e,(t>) (12) 

and 

Rni(r)=N n ,tjt(x ltn r/R) (13) 

In these equations Yi m (9, </>) are the well-known spherical harmonics, ji(z) is a spherical Bessel function and xi iU its 
nth zero[l6(. 

IV. RESULTS AND DISCUSSION 

TableUshows the total, kinetic and potential energies for the ground state for several values of the box radius R. The 
first, second, and third rows for each entry shows perturbation and variational results for the MNC and variational 
results for the CNC, respectively. For simplicity, we restrict our MNC variational calculations to the trial function 

<p(r e , r n ) = (l - ^) (l - ^) e -<*re-Pr„- 7 r ^ + ^ + ^ + ^ (M) 

In the CNC we consider Eq. ((9]) with N = 4 that yields as much as five-digits accuracy when compared with the 
higly accurate results of Aquino et al[l7f. On the other hand, the approximate variational MNC energy is expected to 
be less accurate because the basis set in Eq. (fT4f is comparatively smaller. However, present results are considerably 
more accurate than previous ones for this modelQ, 

First-order perturbation theory is acceptable in the strong-coupling region R < 1 au. We appreciate that the 
perturbation result is systematically greater than the variational one that is known to be a rigorous upper bound. 
This fact clearly shows that the latter is more accurate even for R = 0.1 au where the perturbation result is expected to 
be fairly accurate. TableUshows that perturbation theory underestimates the nuclear kinetic energy and overestimates 
the potential energy while, on the other hand, provides a quite reasonable estimate of the electronic kinetic energy. 

It is also interesting to compare the variational results for the MNC and CNC. The MNC energy is greater than 
the CNC one for all values of R. The reason is that the CNC energy is smallest when the nucleus is located at the 
center of the box that is the particular case chosen here. This effect was discussed earlier by means of a less accurate 
trial function for the MNC 0,0. TableUshows that both the kinetic and potential energies are greater for the MNC. 

Table [TT] shows the variational parameters for the trial functions (fT4|) and ([9|) as well as the expectation values of 
r e , r„ and r. Since the magnitude of the nonlinear variational parameter a is negligible for all values of R we have 
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set it equal to zero. On the other hand, f3 is quite large for small R and decreases as R increases. The remaining 
nonlinear parameter 7 decreases with R reaches a minimum and then increases asymptotically towards the free-atom 
value. Those exponential parameters suggest that the nucleus is localized about the box center whereas the electron 
is localized about the nucleus (though not so strongly because (3 > 7 for all R). The fact that (r n ) < (r e ) ~ R/2 
supports this conjecture. Table ITT1 also shows that (r) MNC ~ ( r )cNC- 
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TABLE I: Total, kinetic and potential energies for the confined hydrogen atom for several values of the box radius. The entries 
stand for MNC perturbation theory, MNC variational method and CNC variational method, respectively 



R 


E 


<T> <T e > <T„> (V) 


0.1 


475.88825 
473.84272 
468.99313 


493.74898 493.48022 0.26876 -17.86073 
497.55784 495.52046 2.03739 -23.71513 
493.59225 - - -24.59911 


0.2 


114.50688 
112.47785 
111.07107 


123.43724 123.37006 0.06719 -8.93037 
124.54287 123.91907 0.62381 -12.06502 
123.48629 - - -12.41522 


0.3 


48.90742 
47.28928 
46.59279 


54.86100 54.83114 0.02986 -5.95358 
55.44914 55.13094 0.31820 -8.15987 
54.94902 - - -8.35623 


0.4 


26.39413 
25.06003 
24.63398 


30.85931 30.84251 0.01680 -4.46518 
31.26099 31.06244 0.19855 -6.20096 
30.96278 - - -6.32881 


0.5 


16.17781 
15.03997 
14.74805 


19.74996 19.73921 0.01075 -3.57215 
20.06358 19.92556 0.13802 -5.02361 
19.86217 - - -5.11412 


0.6 


10.73846 
9.74237 
9.52774 


13.71525 13.70778 0.00747 -2.97679 
13.98064 13.87803 0.10261 -4.23827 
13.83361 - - -4.30586 


0.7 


7.52498 
6.63550 
6.46994 


10.07651 10.07102 0.00548 -2.55153 
10.31303 10.23317 0.07986 -3.67754 
10.19984 - - -3.72990 


0.8 


5.48224 
4.67560 
4.54339 


7.71483 7.71063 0.00420 -2.23259 
7.93312 7.86889 0.06424 -3.25753 
7.84251 - - -3.29912 


0.9 


4.11114 
3.37058 
3.26219 


6.09567 6.09235 0.00332 -1.98453 
6.30211 6.24913 0.5298 -2.93153 
6.22741 - - -2.96522 


1.0 


3.15142 
2.46468 
2.37399 


4.93749 4.93480 0.00269 -1.78607 
5.13617 5.09161 0.04456 -2.67149 
5.07314 - - -2.69915 


1.1 


1.81475 
1.73761 


4.27426 4.23620 0.03807 -2.45952 
4.22006 - - -2.48245 


1.2 


1.33794 
1.26931 


3.64094 3.61044 0.03049 -2.30300 
3.57213 - - -2.30282 


1.3 


0.97551 
0.91704 


3.12685 3.10009 0.02676 -2.15133 
3.06876 - - -2.15172 


1.4 


0.69758 
0.64711 


2.71961 2.69594 0.02367 -2.02203 
2.67018 - - -2.02307 


1.5 







R 


E 


<T) <T e ) <T„> (V) 


2.0 


-0.09946 
-0.12500 


1.43084 1.41801 0.01283 -1.53031 
1.41016 - - -1.53515 


3.0 


-0.41078 
-0.42397 


0.78067 0.77452 0.00615 -1.19145 
0.77206 - - -1.19603 


4.0 


-0.47520 
-0.48327 


0.59178 0.58854 0.00325 -1.06699 
0.58486 - - -1.06813 


5.0 


-0.49138 
-0.49642 


0.53155 0.52961 0.00194 -1.02293 
0.52453 - - -1.02095 


6.0 


-0.49615 
-0.49928 


0.51194 0.51064 0.00130 -1.00809 
0.50635 - - -1.00563 


7.0 


-0.49784 
-0.49986 


0.50515 0.50420 0.00096 -1.00299 
0.50148 0.50148 - -1.00135 


8.0 


-0.49857 
-0.49997 


0.50256 0.50180 0.00075 -1.00113 
0.50032 0.50032 - -1.00029 


9.0 


-0.49894 
-0.50000 


0.50145 0.50082 0.00062 -1.00039 
0.50006 0.50006 - -1.00006 


10.0 


-0.49916 
-0.50000 


0.50070 0.50013 0.00057 -0.99985 
0.50002 0.50002 - -1.00002 



TABLE II: Variational parameters for the trial functions (|14[) and (|9} {N = 4) and some expectation values 



R 


ft 7 Cl C2 C3 C4 

8 do d\ d2 dz d± 


(r e ) (r„) (t)mnc 
(t)cnc 


0.1 


71.386 12.285 0.056 1.000 -0.210 0.446 

4.915 0.001 0.014 0.004 -1.154 2.836 


0.050 0.017 0.052 
0.050 


0.2 


40.494 6.305 0.120 1.000 -0.203 0.514 

1.068 0.095 0.531 -2.000 -1.603 -0.639 


0.099 0.031 0.102 
0.099 


0.3 


29.340 4.287 0.195 1.000 -0.203 0.588 

1.321 0.140 0.551 -0.946 -1.275 -0.567 


0.148 0.044 0.152 
0.147 


0.4 


23.400 3.267 0.286 1.000 -0.207 0.670 

1.267 0.261 0.767 -0.768 -1.127 -0.513 


0.196 0.055 0.200 
0.195 


0.5 


19.647 2.651 0.397 1.000 -0.214 0.767 

1.217 0.419 0.977 -0.593 -0.974 -0.437 


0.243 0.066 0.248 
0.242 


0.6 


17.031 2.237 0.535 1.000 -0.225 0.883 

1.156 0.617 1.171 -0.446 -0.848 -0.333 


0.290 0.077 0.295 
0.289 


0.7 


15.088 1.939 0.694 0.974 -0.233 1.000 

1.095 0.797 1.253 -0.293 -0.706 -0.216 


0.337 0.087 0.341 
0.335 


0.8 


13.578 1.713 0.781 0.825 -0.215 1.000 

1.041 0.958 1.262 -0.163 -0.586 -0.126 


0.382 0.097 0.387 
0.380 


0.9 


12.365 1.536 0.867 0.685 -0.197 1.000 

1.001 1.091 1.225 -0.052 -0.485 -0.065 


0.427 0.107 0.431 
0.425 


1.0 


11.365 1.393 0.954 0.551 -0.182 1.000 

1.011 0.972 0.980 0.073 -0.330 -0.033 


0.471 0.117 0.475 
0.468 


1.1 


17.414 0.480 0.025 -0.091 1.000 0.095 

1.212 0.885 1.000 0.233 -0.074 -0.124 


0.515 0.123 0.518 
0.511 


1.2 


16.215 0.481 0.028 -0.095 1.000 0.099 

1.216 0.948 1.000 0.272 -0.022 -0.099 


0.557 0.132 0.560 
0.553 


1.3 


15.168 0.482 0.031 -0.100 1.000 0.103 

1.224 0.993 0.989 0.298 0.019 -0.076 


0.599 0.141 0.602 
0.595 


1.4 


14.243 0.485 0.033 -0.104 1.000 0.107 

1.264 1.000 0.983 0.328 0.071 -0.061 


0.640 0.150 0.642 
0.635 


1.5 


13.418 0.488 0.037 -0.109 1.000 0.112 

1.279 1.000 0.951 0.335 0.093 -0.041 


0.681 0.159 0.681 
0.675 



R 


P 7 Cl C2 C3 C4 

S do di d'2 c?3 di 


(r e ) (r„) {r)MNC 
{t)cnc 


2.0 


10.303 0.510 0.056 -0.131 1.000 0.133 

1.348 1.000 0.841 0.385 0.062 0.045 


0.869 0.204 0.864 
0.859 


3.0 


4.498 0.565 1.000 -0.459 0.312 0.452 

0.922 1.000 0.243 0.091 -0.012 0.012 


1.178 0.319 1.147 
1.153 


4.0 


2.974 0.642 1.000 -0.324 -0.015 0.319 

0.400 1.000 -0.350 0.086 -0.012 0.001 


1.391 0.449 1.320 
1.342 


5.0 


2.159 0.711 1.000 -0.254 -0.063 0.251 

0.424 1.000 -0.375 0.087 -0.012 0.001 


1.539 0.599 1.408 
1.440 


6.0 


1.643 0.768 1.000 -0.200 -0.090 0.201 

0.465 1.000 -0.367 0.079 -0.010 0.001 


1.657 0.761 1.450 
1.481 


7.0 


1.296 0.813i 1.000 -0.157 -0.109 0.165 

0.510 1.000 -0.346 0.068 -0.008 0.000 


1.770 0.930 1.470 
1.495 


8.0 


1.043 0.851 1.000 -0.122 -0.126 0.139 

0.555 1.000 -0.319 0.057 -0.006 0.000 


1.887 1.106 1.480 
1.499 


9.0 


0.841 0.886 1.000 -0.089 -0.142 0.118 

0.601 1.000 -0.287 0.046 -0.004 0.000 


2.018 1.295 1.486 
1.500 


10.0 


0.747 0.910 1.000 -0.067 -0.151 0.102 

0.935 1.000 0.032 0.009 -0.001 0.000 


2.104 1.412 1.490 
1.500 



